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Abstract 

We calculate spectral functions associated with hadronic current correlation functions for vector 
and pseudoscalar currents at finite temperature. We make use of the Nambu-Jona-Lasinio (NJL) 
model with temperature-dependent coupling constants and temperature-dependent momentum 
cutoff parameters. At low energies, good fits are obtained for the spectral functions that were 
extracted from lattice data by means of the maximum entropy method (MEM). Our model has 
two parameters which are used to fix the magnitude and position of the large peak seen in the 
spectral functions. With those two parameters fixed, we obtain a satisfactory fit to the width of the 
peak. The model then also reproduces the energy of a second peak seen in the spectral functions. 
In the case of the pseudoscalar spectral function, the calculated peak is about 20 percent higher 
than that found for the spectral function obtained from the lattice data. However, it appears 
that the second peak is a lattice artifact [ P. Petreczky, private communication ] and our fit to the 
second peak may not be meaningful. We conclude that the NJL model may have a broader range 
of application than previously considered to be the case, if one allows for significant temperature 
dependence of the parameters of the model, as well as rather large values of the momentum cutoff 
parameter. Our treatment of temperature-dependent coupling constants and cutoff parameters is 
analogous to the procedure introduced by R. Casalbuoni, R. Gatto, G. Nardulli, and M. Ruggieri, 
[Phys. Rev. D 68, 034024 (2003) ], who make use of the NJL model at finite density and find that 
they need to use the density-dependent coupling constants and density-dependent cutoff parameters 
to study matter at high density. 

PACS numbers: 12.39.Fe, 12.38.Aw, 14.65.Bt 
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I. INTRODUCTION 



In a member of recent works [ 1- 3] we have calculated various hadronic correlation func- 
tions and compared our results to results obtained in lattice simulations of QCD [4-6]. The 
lattice results for the correlators, G(t, T), may be used to obtain the corresponding spectral 
functions, a(u,T), by making use of the relation 

PCX) 

G(t,T)= diuap(uj,T)K(T,uj,T) , (1.1) 
Jo 

where 

cosh[c(r-l/2T)] 
K ( T > U > T > = sinh(c/2T) • (L2) 

The procedure to obtain a(ou, T) from the knowledge of G(t, T) makes use of the maximum 
entropy method (MEM) [7-9], since G(t,T) is only known at a limited member of points. 

In our previous work we have made use of the Nambu-Jona-Lasinio (NJL) model. The 
Lagrangian of the generalized NJL model we have used in our studies is 



C = g(* 7 - m°)q + -f ]T[(gA l g) 2 + (qi^X'q) 2 ] (1.3) 

i=0 

-^E^7^) 2 + (^757^) 2 ] 

i=0 

G 

+ ^{det[g(l + X 5 )q] + det[g(l - A 5 )g]} + C Co nf 

Here, m° is a current quark mass matrix, m° = diag{m^,m° d ,m° s ). The \ are the Gell- 
Mann (flavor) matrices and A = a/2/31, with I being the unit matrix. The fourth term is 
the 't Hooft interaction and Cconf represents the model of confinement used in our studies of 
meson properties. For the present work we neglect the 't Hooft interaction and Cconj- Thus, 
there are essentially three parameters to consider, Gs, Gy and a Gaussian cutoff parameter 
a, which restricts the momentum integrals through a factor exp[— k 2 /a 2 }. When we use the 
NJL model to study matter at finite temperature, we introduce the temperature-dependent 
parameters Gs(T), Gy(T) and ct(T). These parameters are adjusted to obtain fits to the 
spectral functions a{u,T) for T/T c = 1.5 and T/T c = 3.0, which are the values of T/T c 
studied in the lattice simulations of QCD that we consider in this work [10]. 
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TABLE I: Parameters used in the calculation of the results reported in Figs. 1 and 2 are given. 



T/T c 





1.5 


3.0 


G s ( GeV" 2 ) 


13.49 a 


4.00 


1.50 


G v ( GeV" 2 ) 


11.46° 


7.50 


2.90 6 


a{ GeV" 2 ) 


0.605 a 


0.850 


1.50 



a) Values used in Ref . [ 20 ] , 

b) For the calculation reported in Fig. 3, this value is changed to Gy = 2.30 GeV -2 . 

Our application of a generalized NJL model for the calculation of temperature-dependent 
hadronic correlation functions has been described in detail in our earlier work [ 1-3 ] . For 
ease of reference, we include the relevant material in the Appendix of this work. 

In Section II we present the values obtained in our analysis of the pseudoscalar and vector 
spectral functions. (The analysis of the axial-vector and scalar spectral functions is quite 
similar to that given in the Appendix for the vector and pseudoscalar spectral functions.) 
In Section III we provide some further discussion of our results and some conclusions. 

II. CALCULATION OF HADRONIC CURRENT SPECTRAL FUNCTIONS 

In this section we present results obtained for temperature-dependent hadronic current 
correlation functions making use of the formalism reviewed in the Appendix. Our results 
may be compared to those presented in Fig. 4 of Ref. [10]. 

Our calculations are made using the parameters given in Table 1. In Fig. 1 we present 
values for the pseudoscalar spectral function divided by the square of the energy, o /uo 2 . The 
solid line in Fig. 1 is obtained when T/T c = 1.5 and the dashed line is for T/T c = 3.0. In 
all our calculations we use a Gaussian cutoff, exp[— A; 2 /a; 2 ], when calculating the vacuum 
polarization functions appearing in the denominator of the spectral function. [ See the 
Appendix. ] The numerator of the spectral function is calculated without a cutoff, so that 
our result goes over to the perturbative result at large energies. Since the denominator of 
the spectral function is calculated using our version of NJL model, a cutoff appears naturally 
in the calculation of the polarization function appearing there. 

We note that, once the two parameters of the model (Gs and a) are fixed, the width of the 
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FIG. 1: Calculated values of ap(T)/uj 2 are shown as a function of lo/T for the model described in 
the Appendix. [See Table 1.] Comparison may be made to Fig. 4 of Ref. [10] or Fig. 4 of the 
present work. Here the solid line corresponds to T/T c = 1.5, while the dashed line is calculated for 
T/T c = 3.0. 

main peak of the pseudoscalar spectral function seen in Fig. 1 is somewhat larger than that 
seen in Fig. 4. The position of the second broad peak seen in the spectral function obtained 
in the lattice study at about 3 GeV for T/T c = 1.5 and at about 6 GeV for T/T c = 3.0 is 
satisfactory. For T/T c = 1.5 the height of the second peak calculated in our work is about 
25 percent larger than that seen in Fig. 4 of Ref. [10]. (We recall that the second peaks 
obtained in the lattice simulations are thought to be unphysical. ) 

Our results for the vector spectral function is given in Fig. 2, where we have use the 
parameters of Table 1. In Fig. 4 of Ref. [10] the height of the peak of the vector spectral 
function for T/T c = 3.0 is significantly lower than the peak for T/T c = 1.5, however, the 
error is quite large in this case so that our result given in Fig. 2 is still consistent with the 
lattice data. In Fig. 3, we use a smaller value of Gy to lower the peak for the T/T c = 3.0 
curve. [See Table 1.] Both the results given in Figs. 2 and 3 are consistent with the values 
extracted from lattice data for T/T c = 3.0, given the large theoretical error in this case. 

As noted in the abstract, the second peak seen in the spectral function is thought to be a 
lattice artifact [ 11 ]. Therefore, the fact that we predict the energy of the second peak may 
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FIG. 2: Calculated values of ay(T)/u> 2 are shown as a function of w/T. [See Table 1.] Here the 
solid line corresponds to T/T c = 1.5, while the dotted line is calculated for T/T c = 3.0. ( See Fig. 4 
of Ref. [ 10 ] or Fig. 5 of the present work ) The theoretical value of ay (T) defined in the appendix 
is here multiplied by 3/4 to correspond to the normalization used in the literature [4-6, 10]. [See 
Eqs. (A30) and (A31)]. 

not be meaningful. 

III. CONCLUSION 

The dynamics of the quark-gluon plasma is a topic of great interest. It has been found 
that the plasma has a number of features that cannot be described in perturbation theory. 
The appearance of resonances in the temperature-dependent spectral functions is one of 
these features. It may be noted in Fig. 5 of Ref. [ 10 ] that resonances appear in the scalar, 
pseudoscalar, vector, and axial-vector spectral functions at about the same energy, which 
indicates that the state-dependent quark interaction plays a less important role in separating 
these states in the deconfined phase than in the confined phase. For example, it is well 
known that the pion loses its character as a Goldstone boson and becomes (approximately) 
degenerate with the sigma meson in the deconfined phase. 

It is of interest to note that the NJL model may be extended to calculate resonances in 
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0.18 -| 
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co/T 

FIG. 3: Calculated values of av(T)/oo 2 are shown. The value of Gy for the dashed line (T/T c = 
3.0 ) is changed from the value of 2.9 GeV" 2 used for Fig. 2 to 2.3 GeV~ 2 . See Fig. 4ofRef. [10] 
or Fig. 5 of the present work. Note the large theoretical error for the vector spectral functions in 
Ref. [ 10 ], which makes the height of the first peak somewhat uncertain. (See the comment on the 
normalization of oy(T) given in the caption of Fig. 2.) 

0.18 -, 
0.16- 
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FIG. 4: The spectral function a/u 2 for pseudoscalar states obtained by MEM [10] is shown. The 
solid line is for T/T c = 1.5 and the dashed line is for T/T c = 3.0. 
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0.18 -, 




to/T 

FIG. 5: The spectral function a/u 2 for vector states obtained by MEM [10] is shown. See the 
caption of Fig. 4. 

the deconfined phase. In order to describe these resonances we have made the parameters 
of the model, Gs, Gy and cv, temperature-dependent. We see that our model provides 
reasonable values for the widths of the resonances and also predicts the position of a second 
broad resonance seen in the spectral functions. (However, the latter resonance is thought to 
be a lattice artifact, as noted earlier.) 

The temperature-dependent coupling constants and cutoff parameters of our work are 
analogous to the corresponding density-dependent parameters introduced in Refs. [ 12 ] 
and [13]. Further study of models with temperature-dependent and density-dependent 
parameters are of interest and a general theoretical formalism for the introduction of such 
dependencies should be considered. 

In future work we also hope to calculate the spatial correlators in the deconfined phase. 
Some data obtained in lattice studies for such correlation functions are given in Ref. [10]. 

APPENDIX 

For ease of reference, we present a discussion of our calculation of hadronic current 
correlators taken from Ref. [ 3 ] . The procedure we adopt is based upon the real-time finite- 
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temperature formalism, in which the imaginary part of the polarization function may be 
calculated. Then, the real part of the function is obtained using a dispersion relation. The 
result we need for this work has been already given in the work of Kobes and Semenoff [14]. 
(In Ref. [ 14] the quark momentum is k and the antiquark momentum isk — P. We will adopt 
that notation in this section for ease of reference to the results presented in Ref. [14].) With 
reference to Eq. (5.4) of Ref. [14], we write the imaginary part of the scalar polarization 
function as 

w « = i( WV^V- (^L m ) (Ai) 

{(1 - rai(fc) - n 2 (k))5(P° - E^k) - E 2 (k)) 
-MA;) - n 2 (k))5(P° + E 1 (k) - E 2 {k)) 
-M*0 - ^(kMP - E^k) + E 2 (k)) 
_(1 _ ni (k) _ n 2 (k))5(P° + E 1 (k) + E 2 (k))} . 

Here, E±(k) = [k 2 + m^T)] 1 / 2 . Relative to Eq. (5.4) of Ref. [14], we have changed the 
sign, removed a factor of g 2 and have included a statistical factor of 2N C , where the fac- 
tor of 2 arises from the flavor trace. In addition, we have included a Gaussian regulator, 
exp[ — k 2 /a 2 ]. The value a = 0.605 GeV was used in our applications of the NJL model in 
the calculation of meson properties at T = 0. We also note that 

"i(*) = e /»frW + l ' ( A2 ) 

and 

H2 ^ = e 0E 2 (k) + i • ( A3 ) 

For the calculation of the imaginary part of the polarization function, we may put k 2 = 
m\(T) and (k — P) 2 = m|(T), since in that calculation the quark and antiquark are on- 
mass-shell. In Eq. (Al) the factor (3 S arises from a trace involving Dirac matrices, such 
that 

(3 S = -Tr[(^ + m 1 )(^-f + m 2 )] (A4) 
= 2P 2 -2(m 1 +m 2 ) 2 , (A5) 
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where m-i and m 2 depend upon temperature. In the frame where P — 0, and in the case 
mi = m 2 , we have (3s = 2-P 2 (1 — 4m 2 /P 2 ). For the scalar case, with m x = m 2 , we find 

52 / A lm 2fr r s\ 3/2 



(A6) 



where 



k 2 = I l- m 2 (T). (A7) 

For pseudoscalar mesons, we replace (3$ by 

P P = -Tr[27 5 (^ + m 1 )2 75 (^-f + m 2 ) ] (A8) 
= 2P 2 -2(mi -m 2 ) 2 , (A9) 

which for m\ = m 2 is (3p = 2P 2 in the frame where P = 0. We find, for the n mesons, 

l mM P\T) = ^ (l - ^2) 1/1 'e-»W [I _ 2„ lW ], ( A10 ) 

where A; 2 = Pq/A — m 2 (T), as above. Thus, we see that the phase space factor has an 
exponent of 1/2 corresponding to a s-wave amplitude. For the scalars, the exponent of the 
phase-space factor is 3/2, as seen in Eq. (A6). 
For a study of vector mesons we consider 

^ = Tr[ 7M (^ + mi ) lu (ft - f + m 2 ) ] , (All) 

and calculate 

=4[P 2 -m 2 -m 2 + 4m 1 m 2 ], (A12) 

which, in the equal- mass case, is equal to 4P 2 + 8m 2 (T), when mi = m 2 and P = 0. This 
result will be needed when we calculate the correlator of vector currents. Note that, for 
the elevated temperatures considered in this work, m u (T) = m d {T) is quite small, so that 
4P 2 + 8m 2 (T) can be approximated by 4P 2 , when we consider the vector current correlation 
functions. In that case, we have 

Im J V (P 2 , T) ~ |lm J P (P 2 , T) . (A13) 
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At this point it is useful to define functions that do not contain that Gaussian regulator: 

Im J P (P 2 , T) = ^2 (l - V2 [ 1 - 2m (k) } , (A14) 

and 

Im J V (P 2 , T) = (i _ ^Tj V2 [ 1 - 2 ni (k) ] , (A15) 

For the functions defined in Eq. (A14) and (A15) we need to use a twice-subtracted dispersion 
relation to obtain Re J P (P 2 ,T), or Re Jy(P 2 ,T). For example, 

P 2 

Re J P (P 2 , T) = Re J P (0, T) + [ Re J P (P 2 , T) - Re J P (0, T) ] (A16) 
, P\P 2 ~ Pj) ds _ lraJ P (s,T) 

71 A 



W) s(P 2 - s)(P 2 - s) ' 

where A 2 can be quite large, since the integral over the imaginary part of the polarization 
function is now convergent. We may introduce J P (P 2 ,T) and Jy(P 2 ,T) as complex func- 
tions, since we now have both the real and imaginary parts of these functions. We note that 
the construction of either Re J P (P 2 ,T), or Re Jy(P 2 ,T), by means of a dispersion relation 
does not require a subtraction. We use these functions to define the complex functions 
J P (P 2 ,T) and J V (P 2 ,T). 

In order to make use of Eq. (A16), we need to specify J P (0) and Jp(P 2 ). We found it 
useful to take P 2 = -1.0 GeV 2 and to put J P (0) = J P (0) and Jp(P 2 ) = Jp(P 2 )- The 
quantities Jy(0) and Jy(P 2 ) are determined in an analogous function. This procedure in 
which we fix the behavior of a function such as ReJy(P 2 ) or ReJy(P 2 ) is quite analogous 
to the procedure used in Ref. [15]. In that work we made use of dispersion relations to 
construct a continuous vector-isovector current correlation function which had the correct 
perturbative behavior for large P 2 — > — oo and also described the low-energy resonance 
present in the correlator due to the excitation of the p meson. In Ref. [ 15 ] the NJL model 
was shown to provide a quite satisfactory description of the low-energy resonant behavior 
of the vector-isovector correlation function. 

We now consider the calculation of temperature-dependent hadronic current correlation 
functions. The general form of the correlator is a transform of a time-ordered product of 
currents, 

tC(P 2 ,T) = J d A xe iP x « T(j(x)j(0)) » , (A17) 
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where the double bracket is a reminder that we are considering the finite temperature case. 

For the study of pseudoscalar states, we may consider currents of the form jp,i(x) = 
q(x)i'j 5 \ z q(x), where, in the case of the n mesons, i — 1,2 and 3. For the study of scalar- 
isoscalar mesons, we introduce js,i( x ) — q{x)\ l q(x), where % — for the flavor-singlet current 
and % — 8 for the flavor-octet current [16]. 

In the case of the pseudoscalar-isovector mesons, the correlator may be expressed in terms 
of the basic vacuum polarization function of the NJL model, Jp(P 2 , T) [ 17-19]. Thus, 

c ^^^ ^Wii ' (A18) 

where Gp(T) is the coupling constant appropriate for our study of ir mesons. We have 
found Gp(T) = 13.49 GeV~ 2 by fitting the pion mass in a calculation made at T = 0, with 
m u = irid = 0.364 GeV. The result given in Eq. (A18) is only expected to be useful for small 
P 2 , since the Gaussian regulator strongly modifies the large P 2 behavior. Therefore, we 
suggest that the following form is useful, if we are to consider the larger values of P 2 . 



C P (P 2 ,T) 
P 2 



Jp(P 2 ,T) 
P 2 



' (A19) 



l-Gp(T)Jp(P*,T)' 



(As usual, we put P = 0.) This form has two important features. At large Pq, 
Im Cp(P , T)/P 2 is a constant, since Im J p (Pq,T) is proportional to P 2 . Further, the 
denominator of Eq. (A19) goes to 1 for large P 2 . On the other hand, at small P 2 , the 
denominator is capable of describing resonant enhancement of the correlation function. As 
we have seen, the results obtained when Eq. (A19) is used appear quite satisfactory. ( We 
may again refer to Ref. [15], in which a similar approximation is described.) 

For a study of the vector-isovector correlators, we introduce conserved vector currents 
j[t,i( x ) — Q( x )l 'n\q( x ) with i=l, 2 and 3. In this case we define 

jr(P 2 ,T) =(g^- Jv{P\T) (A20) 

and 

CT(P 2 ,T) =[g^- -p^J CV(P 2 ,T) , (A21) 

taking into account the fact that the current j^i(x) is conserved. We may then use the fact 
that 

Jv(P 2 ,T)= 1 -g, I/ J^(P 2 ,T) (A22) 
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and 



Im J V (P 2 ,T) 



Pg + 2ml(T) 
An 



-ImJ P (P 2 ,T). 



1 - 



p2 



1/2 



(See Eq. (A7) for the specification of k — \k\) We then have 

1 



C V (P 2 ,T) = J V (P 2 ,T) 



1-G V (T)J V (P*,T)' 



e -k 2 /* 2 [ 1 _2n 1 (k)} (A23) 

(A24) 



(A25) 



where we have introduced 



Im J V (P 2 ,T) 



n 2 + 2m^(T) 
An 



±ImJ P (P 2 ,T). 



4m 2 (T) 



p2 



1/2 



l-2 ni (k) 



In the literature, tu is used instead of P [ 4-6 ] . We may define the spectral functions 



and 



ay (u, T) = — Im Cv(uj, T) , 
n 



<tp(u,T) = -1mC P (u,T) , 
n 



(A26) 
(A27) 

(A28) 

(A29) 



Since different conventions are used in the literature [4-6], we may use the notation 
ffp(w, T) and a v (uj, T) for the spectral functions given there. We have the following relations: 



and 



a P (u,T) =a P (u,T), 



a v (uj,T) 3 

2 = i°v(u,T), 



(A30) 



(A31) 



where the factor 3/4 arises because, in Refs. [4-6], there is a division by 4, while we have 
divided by 3, as in Eq. (A22). 
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